Microfluidic velocimetry reveals spatial cooperativity in soft glassy flows 
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Amorphous glassy materials of diverse nature 
— concentrated emulsions, granular materials, 
pastes, molecular glasses display complex flow 
properties, intermediate between solid and liquid, 
which are at the root of their use in many appli- 
cations [l\ [2l |3]. A classical feature, well docu- 
mented yet not really understood, is the very non- 
linear nature of the flow rule relating stresses and 
strain rates [4l IB]. Using a microfluidic velocime- 
try technique, we characterize the flow of thin 
layers of concentrated emulsions, confined in gaps 
of different thicknesses by surfaces of different 
roughness. Beyond the classical non-linearities of 
the rheological behaviour, we evidence finite size 
effects in the flow behaviour and the absence of an 
intrinsic local flow rule. In contrast, a rather sim- 
ple non-local flow rule is shown to account for all 
the velocity profiles. This non-locality of the dy- 
namics is quantified by a length, characteristic of 
the cooperativity of the flow at these scales, that 
is unobservable in the liquid state (lower concen- 
trations) and that increases with concentration 
in the jammed state. Beyond its practical impor- 
tance for applications involving thin layers, e.g. 
coatings, our assessment of non-locality and co- 
operativity echoes observations on other glassy, 
jammed and granular systems, suggesting a pos- 
sible fundamental universality. 



Glasses and jammed matter are at the heart of many 
challenging questions in condensed matter physics. The 
glass transition and the physical and mechanical prop- 
erties of materials in the glassy or jammed state, from 
granular materials to dense colloidal suspensions, remain 
puzzle to a large extent H H |5J El |S]. The 
transition itself is characterized by a dramatic slowing 
down of the system, associated with a dynamical ar- 
rest of the microscopic structure [3]. While a large 
amount of work has been dedicated in recent years 
to demonstrating and characterizing dynamical hetero- 
geneities in granular, glassy and glass forming materials 
[lEnilfflEIimUlEnEE], the question of whether 



and how these dynamical heterogeneities influence the 
flow behaviour remains [17l [TSJ US] . Indeed one expects 
that such cooperative effects, which are in essence dy- 
namical processes, should manifest themselves in the flow 
dynamics of jammed materials. This question is essen- 
tial from the perspective of both fundamental science and 
applications: flows of jammed and glassy materials and 
glassy films are ubiquituous in nature and industry, e.g. 
in granular flows, coatings with thin polymer films, lubri- 
fying process in solid friction, in the food and cosmetic 
industries, and even in pedestrian dynamics |20j . 

In the present work, we aim to develop a constitu- 
tive law governing flows of jammed materials. We take 
advantage of a local velocity measurement technique to 
follow the local flow behaviour of a film of a confined 
soft glassy material (here a concentrated emulsion). We 
demonstrate the existence of finite-size effects in the flow 
dynamics that cannot be apprehended within our current 
understanding of the behaviour of yield stress fluids. 

As schematized in Fig. 1, we probe the flow of 
a jammed emulsion in various geometries and confine- 
ments: shear or pressure driven planar flow; wide gap 
Couette cell (one centimer gap) or narrow microchan- 
nel (from tens to hundreds of microns in width). The 
emulsion is made of silicone droplets (« 6.5/im in diam- 
eter) in an index-matched glycerine- water mixture (Sup- 
plementary Information, Method). We are specifically 
interested in the local flow curves of the material, which 
relate the local shear stress a to the local shear rate 7, 
Fig. 1C-D. Such curves can be constructed from the ve- 
locity profiles measured in the jammed emulsion in both 
geometries using local velocimetry techniques (Supple- 
mentary Information, Method), see insets of Fig. 1C-D. 
Indeed in both geometries the stress distribution is known 
from mechanical equilibrium: cr(R) = 2 itHR 2 a ^ a P os i~ 
tion R in the Couette cell (of height H) under a torque T; 
a(z) = ^£-{z — y), at a position z — If from the center- 
line of the microchannel, under an applied pressure drop 
AP between the two ends of the channel of length L [21J. 
Local shear rates 7 are deduced directly from the veloc- 
ity profiles. As shown in Figs.l.C and l.D, the local flow 
curves exhibit very different behaviours depending on the 
confinement of the material. In the wide gap Couette cell, 
all local flow curves obtained for different torques show 
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a perfect superposition (Fig l.C). The flow curve is fur- 
thermore well described by a classical Herschel-Bulkley 
(HB) model, a = a a + A'j 1 / 2 , with cr a 'dynamical' yield 
stress and A a material constant. On the other hand, 
in the narrow microchannels, data for different pressure 
drops AP are scattered in the whole figure and do not 
collapse on a single rheological curve (Fig l.D), see also 
Supplementary Figure SI-4. Accordingly, while a single 
HB model accounts perfectly for all the velocity profiles 
in the wide gap Couette cell (inset of Fig. l.C), such is 
not the case for the narrow microchannels (Supplemen- 
tary Figure SI-8). Furthermore, below the jamming point, 
(f> < <fi c , the local HB rheological model is fully able 
to reproduce the flow behaviour, even in narrow chan- 
nels (Supplementary Figure SI-3). In summary, a first 
key feature emerges from the experimental data: in the 
jammed state, <f> > (f> c , there is no universal local rela- 
tionship between stress and shear rate. 

As a further confirmation of this point, the overall 
shape of the velocity profiles, for the same wall shear 
stress, is shown to vary with the thickness of the chan- 
nel (Fig. [2]) , demonstrating finite-size effects in the flow 
properties of the jammed emulsion. Rougher surfaces are 
shown to induce higher shear rates for the same confine- 
ment and shear stress at the wall. We did furthermore 
verify that the above results do not stem from a struc- 
tural change of the emulsion itself: we checked that nei- 
ther the size distribution of the droplets, nor their shape, 
was affected by the flow and in all the geometries no vari- 
ation of local droplet concentration was measured within 
the gap (Supplementary Method). Altogether, our re- 
sults suggest that the flow does not couple to a structural 
order parameter, so that non-local effects are 'intrinsic' 
to the rheology of the jammed emulsion, in line with sim- 
ilar observations in numerical simulations for glassy flows 

A second crucial feature is the experimental proof of 
a specific surface rheology at the channel walls, distinct 
from the bulk flow behaviour. First, as commonly ob- 
served in flows of concentrated systems, slip effects show 
up at the confining walls [22]. We have characterized 
these effects by measuring the velocity slip-surface stress 
relationship (Fig. |3i), for two different roughness (Sup- 
plementary Method and Supplementary Figure SI-2). 
Strikingly, a remarkable collapse of the data is exhibited 
for various confinements and pressure drops, suggesting 
that this relationship does only depend on the surface 
roughness, see also Supplementary Figure SI-7a. How- 
ever, a careful inspection of Fig. |2]shows that beyond slip 
effects, the shape itself of the flow profiles does depend on 
the surface characteristics of the confining wall, suggest- 
ing that the wall roughness makes the system more fluid. 
In an attempt to rationalize this behaviour we plot the 
wall shear rate as a function of the wall shear stress in 
Fig. [3|d (see also Supplementary Figure SI- 7b). Again a 
striking collapse of all data for various confinements and 
pressure drops is exhibited. We emphasize that such a 
behaviour is not expected a priori since the same wall 



stress can be achieved with different confinements and 
pressure drops. 

On the basis of these experimental facts, we now pro- 
pose a theoretical framework able to rationalize these 
finite-size effects in the rheological behaviour. In concen- 
trated emulsions, and more generally in soft glasses, flow 
occurs via a succession of reversible elastic deformations 
and local irreversible plastic rearrangements associated 
with a microscopic yield stress [2U These 
localized plastic events induce a non-local, long ranged, 
elastic relaxation of the stress over the system. The num- 
ber of plastic rearrangements per unit time, v, plays the 
role of an inverse relaxation time and therefore controls 
the flow and relaxation of the material: a higher rate is 
associated with a more fluid system. Up to an elastic 
modulus G, the rate of rearrangements v is equivalent to 
a "fluidity", /, defined as a = jj, a quantity that has 
been used in recent studies to characterize the rheology 
of yield stress fluids [25]. In the absence of any non- 
local effect, the fluidity would reduce to its bulk value, 
/bulk =7/ Cbuik [note that for the HB expression reported 
above, /bulk oc 7 in the quasistatic limit, as a — » <jq]. 
Now, due to the non-local elastic relaxation of plastic re- 
arrangements ('dynamical cooperativitv')[2~il I2S], a dy- 
namically active region will induce agitation of its neigh- 
bours, and thus a higher rate of plastic rearrangements 
(and vice- versa for a dynamically quiescent region). This 
suggests that non-local effects will occur for the rear- 
rangements rates and consequently for the fluidity, /. In 
order to capture this physical picture, we assume that 
the local fluidity, f(z), obeys a non-local equation in the 
steady state: 



f{z) = /bulk + £ 



2 d 2 f(z) 
d 2 z 



(1) 



In this equation, the non-local term describes how the 
plastic activity spreads spatially over the system due to 
non-local elastic relaxation: it accounts for the coopera- 
tivity involved in the plastic events occuring during flow 
and £ denotes a bulk 'flow cooperativity length'. 

Together with the definition f(z) = j(z)/a(z), Eq. 
([lj constitutes our non-local rheological model for the 
emulsion. An analytical solution of Eq. can be ob- 
tained, giving the fluidity at some position z in terms 
of an integral of the shear-stress over the whole channel 
width. For convenience we solved Eq. numerically 
with f(z) = 7waii/cwaii as a given boundary condition 
at the walls and / = in the centerline. Velocity pro- 
files are deduced by integrating the computed local shear 
rate j(z). Let us first note that whatever the boundary 
condition at the wall the solution of Eq. (TXT> can not re- 
duce to its local form, fbuik{z) = fbuik(v(z)), due to the 
diffusive, non-local term. 

The only remaining parameter is the bulk 'flow coop- 
erativity length' £ in Eq. (TXT) . A crucial and surpris- 
ing finding is that a unique (constant) length £ accounts 
for all experimental data for the flow profiles and local 
flow curves for a given emulsion, Figs. l.D and 2, inde- 
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pendently of the pressure drop, confinement and surface 
nature (rough or smooth) . That so much data can be fit- 
ted using a single and constant length £ is quite remark- 
able. Typical comparisons are shown in Figs. [Tp and [2] 
for both the velocity profiles and the corresponding local 
flow curves (sec also Supplementary Figures SI-5, SI-6). 
We mention that alternative non-local rheological mod- 
els, involving e.g. non-local terms in the flow curve |27] 
are unable to reproduce all measured flow profiles, Sup- 
plementary Figure SI-9. 

Futhermore, a similar, very good agreement has been 
obtained for all (jammed) emulsions investigated, with 
various volume fractions and polydispersities (Supple- 
mentary Figures SI-5 and SI-6). Fig. [4] reports the vari- 
ation of the measured flow cooperativity length £ as a 
function of the volume fraction for emulsions with two 
different dispersities. Below the jamming concentration 
< 0c we could not measure finite-size effects in the flow 
profiles within the experimental uncertainty (Supplemen- 
tary Figure SI-3), so that £ » in this case. A strong 
increase of £ is observed above the jamming threshold, in 
line with the appearance of elastic and yield properties 
of the materials in this regime. 

Non-local effects in the flow properties may be a pri- 
ori expected at the microscopic scale [27] . However, we 
measure a 'flow cooperativity length" only in the jammed 
state, providing a direct and novel connection between 
non-local effects and the jamming transition. Further- 
more, we emphasize that this length is associated with 
correlations of a dynamical quantity (the fluidity, i.e. re- 



arrangements) and not with static structural order. Typ- 
ically for the volume fraction range studied here, £ varies 
between and 30 microns, to compare with the droplet 
size d « 6fim, giving roughly a cooperatively rearrang- 
ing region involving up to a few hundreds of particles. 
These values are close to estimates for the dynamical 
heterogeneities measured in glassy materials |101 111] or 
in jammed granular materials |14l [TB . 28J. Furthermore, 
similar observations of non-locality have been reported in 
granular flows close to the jamming transition [HI [29] . 
suggesting further universal characteristics. 

However, we emphasize that the present length char- 
acterizing flow cooperativity differs fundamentally from 
that which characterizes dynamical heterogeneities in- 
volved in spontaneous fluctuations. Indeed, while these 
dynamical heterogeneities have been measured to ex- 
hibit a maximum at the jamming point |10[ 130] - the 
present flow cooperativity length £ is non- vanishing in the 
jammed region only, and increases when going deeper in 
the jammed phase. This suggests that the flow behaviour 
of glassy and jammed materials involves physical mecha- 
nisms distinct from those associated with the glass tran- 
sition. We expect that the flow cooperativity length £ 
measures the zone of influence of localized plastic events 
occuring during flow and not the size of mobile regions. 
Our results thus provides a novel perspective for devel- 
oping a theoretical framework to describe the non-local 
flow behaviour of glassy and jammed systems, such as 
soft glasses and granular materials. 
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FIG. 1: Local flow curves measured by IRM and Microfluidic velocimetry. 

Left: Schematic of the experimental set up. (A) wide gap Couette cell (inner and outer cylinder radii, Ri — 4.1cm, and 
R e = 5.9cm respectively, height h — 11cm). (B) Flow in a microchannel and particle imaging velocimetry. PIV : (DM): 
Dichroic mirror; UV : Ultra- Violet lamp; Obj : microscope objective. Flow direction is along x. Right: Local flow curves 
extracted from the measurements of the velocity profiles in both geometries for an emulsion with volume fraction 4> — 0.75 and 
20% polydispersity. The shear rates are obtained by finite differencing of the velocity profiles given in the insets. (C) Local 
flow curves, a(z) versus 7(2), extracted from the velocity profiles measured within a wide gap Couette geometry (see inset). 
The solid line corresponds to the Herschel-Bulkley model with <xo = 11-6 Pa, A = 11.2 Pa.s 1//2 . Inset: Dimensionless velocity 



profiles 



ng) 

V(Ri) 



as a function of the radial coordinate R. From bottom to top the rotation velocity is 5, 10, 20, 50 and 100 rpm. 



(D) Local flow curves, cr(z) versus 7(2), extracted from the velocity profiles measured in a w = 250/im thick microchannel 
with rough surfaces, for various pressure drops as a function of the reduced coordinate z/w (see inset). No overlap of the local 
flow curves is observed. Dashed lines are predictions for the local flow curves at the given AP, as obtained from the non-local 
rheological model of Eq. |l| with a cooperativity length £ = 22.3/xm. The solid line corresponds to the Herschel-Bulkley model 
with o"o = 11.6 Pa, A — 11.2 Pa.s 1 ' 2 . Inset: Corresponding velocity profiles measured for AP equal to 300, 450, 600, 750, 900 
mbar. The length of the channel is L — 14 cm. Solid lines are the velocity profiles predicted by the non-local rheological model 
of Eq. |TJ. Slippage, which is found to occur at the surfaces as detailed in Fig. 3, does not affect the resulting local flow curves 
as plotted in the main figure. 
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FIG. 2: Velocity profiles in microchannels with different gaps w and confining wall roughness, corrected for 
the slip velocity: (a) o rough microchannel with w = 250 /im; (b) x rough microchannel with w = 125 /zm , (c) ■ smooth 
microchannel with w — 112 fim; (d) A smooth microchannel with hi = 56 /im. Solid lines are velocity profiles predicted by 
the non local model of Eq. with a length £ = 22.3/zm. The various experimental velocity profiles correspond to different 
pressure drops AP, tuned to get the same range of wall shear stress <r W aii in the various geometries. From bottom to top the 
values of cr wall are: (a) 27, 41, 55, 68, 82 Pa; (b) 30, 45, 60, 75, 90 Pa; (c) 45, 60, 75, 91 Pa; (d) 48, 65, 82, 97 Pa. The volume 
fraction is <f> = 75%, with 20% polydispersity. Errobars denote s.d. (n=450). 



7 



V (mm.s 1 ) v M (s _1 ) 

s 1 wall 




50 100 50 100 

° wall (Pa) a_ M (Pa) 



wall 




JJLP 



FIG. 3: Rheology at surfaces: (a) Slip velocity at the wall 
as a function of the wall shear stress for various microchannels; 
(b) Shear rate at the confining wall 7 wa ii as a function of 
the wall shear stress own for various microchannels. The 
dashed line is the bulk flow curve (Herschel-Bulkley model 
with o- = 11.6 Pa, A = 11.2 Pa.s 1/2 ). Data gather different 
pressure drops AP and confinements w. Symbols are identical 
to those in Figure [2] Open (filled) symbols correspond to 
rough (smooth) confining surfaces. The volume fraction is 
4> = 75%, with 20% polydispersity. Errobars on the velocity 
denote s.d. (n=450). 



FIG. 4: Flow cooperativity length £ as a function of 
emulsion volume fraction. Open and filled symbols cor- 
responds to emulsions with different polydispersities of 20% 
(open) and 36% (filled) respectively. The corresponding jam- 
ming volume fraction, as determined by the occurence of a 
non-vanishing yield stress, are (f> c = 0.64 (open) and <f> c = 0.68 
(filled). Errobars denote the dispersion of the experimental 
values over various experiments. 



